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Abstract 

We employ Monte Carlo simulations in order to investigate critical behavior of a geometri- 
cally frustrated spin-1 Ising antiferromagnet on a triangular lattice in the presence of a single- ion 
anisotropy. It has been previously found that long-range order can exist in the isotropic system 
with a spin larger than some critical value estimated as 11/2. We show that the presence of the 
single-ion anisotropy can lead to a partial long-range order in the low-temperature region even 
below this critical value, namely for the spin 1, within a certain range of the anisotropy strength. 
At higher temperatures we identify another phase of the Berezinsky-Kosterlitz-Thouless type and 
using a finite-size scaling analysis evaluate the correlation decay exponent. We also study densities 
of various local spin patterns in the respective phases. 
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I. INTRODUCTION 



The triangular lattice Ising antiferromagnet (TLIA) with the spin 1/2 is well known to 
display no long-range ordering down to zero temperature [l| where the spin-correlation func- 
tion decays as power-law [2]. On the other hand, in the ground state a long-range order with 
the partially ordered structure, characterized by two sublattices of opposite magnetizations 
and one sublattice of zero magnetization, can occur if the spin is larger than some critical 
value 3-l5|. The upper bound of this critical value was estimated by the use of Peierls' argu- 
ment 3| as 62 and the precise value was established by Monte Carlo simulations as 11/2. 
The emergence of the long-range order was argued to arise as a result of a drastic change in 
the ground-state degeneracy as well as some other physical quantities of this frustrated spin 
system with the increasing spin magnitude. Nevertheless, t 
by some perturbations, such as an external magnetic field 
which have been shown to lead to long-range ordering even in the highly frustrated spin- 
1/2 system. For Ising models with spin larger that 1/2 a single-ion anisotropy is another 
parameter that may play a crucial role in their critical properties (see, e.g., Ill-ll3|). 

In the present paper we show that this is also the case for the current model and demon- 
strate how the inclusion of the single-ion anisotropy can change the above scenario. We show 
that even a small amount of the single-ion anisotropy leads to a partially ordered phase at 
low temperatures, while different critical phase of the Berezinsky-Kosterlitz-Thouless (BKT) 
type [l4|, characterized by a power-law decay of the spin correlation function, persists at 
higher temperatures. 



le large degeneracy can be hfted 

in nn 

6|-|8| or selective dilution [9|, I10| . 



II. SIMULATION DETAILS 



We consider the Ising model described by the Hamiltonian 

H = -JJ2s^SJ-DJ2sl (1) 

where S'j = ±1, is an Ising spin on the ith lattice site, {i,j) denotes the sum over nearest 
neighbors, J < is the antiferromagnetic exchange interaction parameter, and D is the 
single-ion anisotropy parameter. 

We perform standard Monte Carlo (MC) simulations following the Metropolis dynamics. 
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We simulate the triangular lattice spin systems of the linear size L and apply the periodic 
boundary conditions. For thermal averaging we typically consider = 2 x 10^ MCS (Monte 
Carlo sweeps or steps per spin) after discarding another A^^o = 0.2 x MCS for thermaliza- 
tion. In order to obtain dependencies of various thermodynamic quantities on the reduced 
temperature t = ksT /\J\, we use L = 48 and start simulations from high temperatures with 
random initialization. Then the temperature is gradually lowered with the step At = 0.02 
and the simulations start from the final configuration obtained at the previous temperature. 
In order to estimate the error bars, we run five simulations with different initializations. For 
finite-size scaling analysis (FSS), employed to identify different phases, we use the lattice 
sizes L = 24, 48, 72, 96 and 120. 

We calculate the following quantities. The sublattice magnetization of the sublattice X 
(X=A,B and C) is obtained as Mx = Xliex ^i- Then, for an antiferro magnet, as an order 
parameter it is useful to define the staggered magnetization per site as 



m. 



= {Ms)/L^ = 3(max(MA, Mb, Mc) - min(MA, Mb, Mc))/2L^ (2) 



where (■ ■ ■ ) denotes thermal average. Further, from fluctuations of the above quantities we 
define the specific heat per site c as 



and the staggered susceptibility per site Xs as 



(3) 



In the ground state the spin-correlation function of the TLIA model decays as a power 
law Q 

{SiSj) oc r,7, (5) 

where the exponent rj has been shown to decrease with the spin value from rj = 1/2 for 
spin-1/2 to zero for spin larger than 11/2, for which antiferromagnetic long-range order 
occurs jsl, 01. The value of rj can be estimated by FSS of the order parameter m^, which 
scales as [l5| 



m 



,(L)ocL-''/2 (6) 



or the quantity 
Y = 



m 



which scales as 



Y{L) oc L^-"^. 



i&A jGB k£C 



(7) 



(8) 



III. RESULTS AND DISCUSSION 



Let us first examine the ground-state properties for different values of D. For D = 0, 
the energetic arguments dictate that the ground-state configuration is such that the spins 
on each elementary triangular plaquette sum to ±1. All such plaquettes are energetically 
equivalent with the energy per spin e = J and there is no long-range order among them. 
It is easy to verify that for > the preferred configurations are such that the spins on 
each elementary triangular plaquette again sum to ±1 but zero spin states are not involved. 
Then the energy per spin is e = J — D and the system behaves like a spin-1/2 Ising model 
with no long-range order If D < the configurations with the spin Si = ±1 states on 
one sublattice are suppressed and the partially ordered phase with the antiferromagnetic 
ordering on two sublattices and non-magnetic states of the third sublattice, i.e., the Si = 
states, can occur. Such configurations are characterized by the energy per spin e = J—2D/3. 
However, for D < 3 J/2 the energy becomes positive and therefore the non-magnetic state 
with all spins taking zero value is the ground state. Therefore, the partial long-range order 
can only be expected within 3 J/2 < D < 0. 

In order to confront the behavior of the system with no single-ion anisotropy, which is 
expected to show only the BKT phase transition j^, with the long-range order behavior 
predicted to appear within the range 3J/2 < D < 0, in Fig. [1] we plot the temperature 
dependencies of some relevant thermodynamic quantities for two values of -D/|J| = and 
—1, representing the two cases. The behavior for D/\J\ = is clearly different from that 
observed for D/\J\ = —1. While the former shows only one anomaly at higher temperatures, 
the latter displays an additional anomaly at lower temperatures. In particular, as observed 
in the behaviors of the sublattice magnetizations and the staggered magnetization (order 



parameter), presented in Figs. 1(a) and 1(b), respectively, in both cases some degree of 
ordering is initiated in two sublattices at higher temperatures. However, unlike the case 



with no anisotropy, for D/\J\ = —1 this initial increase in the quantities ttia, m-B and rus 
is followed by another increase to the fully saturated values at lower temperatures. As the 
temperature is lowered, one and two anomalies are also observed in the behavior of the 



internal energy, tending to the expected ground-state values of e/\J\ 
the cases of D/\J\ = and D/\J\ - 



-1 and —1/3, for 



— 1, respectively (Fig. 1(c)). The anomalies in the 
staggered magnetization and the internal energy are reflected in the corresponding number 



of peaks in the staggered susceptibility and the specific heat, shown in Figs. 1(d) and 1(e) 
respectively. 

The above anomalies in the thermodynamic quantities suggest occurrence of transitions 
between different phases. To identify the respective phases, we employ FSS analyses of the 
order parameter and the quantity Y, defined by the relations and ([HD, respectively. 
Temperature variation of the scaling exponent r] is presented in Fig. |2j The value of rj = 0, 
observed for D/\J\ = —1, means that the^sj^stem displays at sufficiently low temperatures 



(below kBT/\J\ ^ 0.35) long-range order [31,1^. On the other hand, for D/\J\ = 0, the value 
of rj remains finite and in the low temperature region levels off at the value of ~ 0.317, in 



agreement with the previous ground-state investigations [3|, |5|. 

Finally, we found it interesting to look into the system's behavior also at microscopic level. 
More specifically, we studied the snapshots of spin configurations at different temperatures 
both visually and quantitatively. For the latter we introduced the quantity, which we termed 
density of local patterns n,, which is defined as a relative number of different spin patterns 
Pi = {Sj, Sk, Si), i = 1 — 14 [l6|, on the elementary triangular plaquettes observed in the 
simulated states. The temperature variations of for different values of the single-ion 



anisotropy strength are shown in Figs. 3(a) 3(c), with all the patterns listed in Fig. 3(a^ 



The snapshots of the respective spin configurations in the ground state are also presented 



in Figs. 3(d) 3(f) [iTj]. For comparison we also included the case with positive value of 
D/\J\ = 1. For D/\J\ = 1, the degenerate spin patterns which sum to ±1 and include 
no zeros, i.e., {±p3, ±p7, ipg}, are equally represented and, in accordance with the above 
ground-state arguments, for t — )• are the only ones present in the system. On the other 
hand, for D/\J\ = —1, due to broken symmetry the system in the partially ordered phase 
chooses one the two configurations involving either the set of patterns {pq, —ps, —Pu} or 
{—^65^85^12}- Again, in agreement with the ground-state arguments, for t — )• the patterns 
of either set are equally represented and there are no other patterns in the system. In 
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between the two phases for D /\J\ = 0, the patterns {±^3, ±P7, ±^9} and {ipe, TPs, TP12} 
mix together since they contribute with the same energy. However, as evidenced from 



both Figs. 3(b) and 3(e), they do not appear with the same probabihties. Namely, the 
patterns {±^3, ±^7, ipg} systematically prevail in all the simulations started from various 
initial states and in the ground state they are almost 2.7-times more represented than the 
patterns {±pq,tP8,TPi2}- This can be understood by looking at larger patterns, namely, 
hexagonal ones formed by a central spin and its six nearest neighbors. Then, only in the 
patterns with the so called "free" central spins (i.e., the spins with the nearest neighbors of 
alternating signs) the state of the central spins is irrelevant with the respect to the energy 
and, therefore, the states +1, —1 and are equally probable. Therefore, the number of the 
triangular patterns {Tpe, TPs, -FP12} is proportional to the number of the hexagonal patterns 



with the "free" central spins in the configurations such as the one shown in Fig. 3(d) , noting 
that the zero-state central spin generates six triangles with the patterns {Tpe, TPs, ^^12}- 



IV. CONCLUSIONS 



We studied the geometrically frustrated spin-1 triangular Ising antiferromagnet in the 
presence of a single-ion anisotropy. We showed that within the range of negative values of the 
anisotropy parameter D/\J\ G (—1.5, 0) the system displays a partial long-range order in the 
low-temperature region and another phase with algebraically decaying correlation function 
at higher temperatures. For D/\J\ < —1.5 the system is in a non- magnetic state and for 
D/\J\ > it behaves like a spin-1/2 model, with no long-range ordering due to frustration. 
The ground-state phase transition at -D/| J| = is characterized by the coexistence of the 
energetically equivalent patterns of the types (±1 ± 1 =F 1) and (0 ± 1 =F 1) at a constant ratio 
of approximately 27 : 10. 
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FIG. 1: Temperature variation of (a) the sublattice magnetizations mA, mB, and mc, (b) the 
staggered magnetization m^, (c) the internal energy e/|J|, (d) the staggered susceptibility Xs and 
(e) the specific heat c, for D/\J\ = —1 and 0. 
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FIG. 2: Temperature variation of the exponent rj, for D/\J\ = —1 and 0. The symbols and the 
error bars respectively represent the mean and the extreme values obtained from Eqs. ([6]) and ([8]). 
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(a)L>/| J| = 1 (b)L'/|J|=0 {c)D/\J\ = -l 




FIG. 3: (Color online) Temperature variation of the densities of local patterns pi = {Sj, S^, Si), 
i = 1 — 14, (upper row) and the corresponding ground-state spin configuration snapshots (lower 
row), for D/\J\ = 1,0 and —1. The red (light gray), blue (dark gray) and white circles represent 
the spin states —1, +1 and 0, respectively. 
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